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Abstract 

The Full (non-truncated) Israel-Stewart theory of bulk viscosity is applied 

to dissipative FRW spacetimes. Dimensionless variables and dimensionless 

equations of state are used to write the Einstein-thermodynamic equations 

as a plane autonomous system and the qualitative behaviour of this system 

is determined. Entropy production in these models is also discussed. 
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I. INTRODUCTION 



In a recent paper JI|, isotropic and spatially homogeneous viscous fluid cosmological 
models were investigated using the Truncated Israel-Stewart |2|-Q theory of irreversible 
thermodynamics to model the bulk viscous pressure. Although it provides a causal and 
stable second order relativistic theory of thermodynamics, the Truncated version of the 
theory can give rise to very different behaviour than the Full Israel-Stewart theory |5|-^J. 
It can be argued that the Truncated theory agrees with the Full theory if one uses, instead 
of the local equilibrium variables, a generalised temperature and thermodynamic pressure 
HIJ. However, there are difficulties in modelling these generalised variables in cosmology. 
Therefore the analysis of [l[ can only be regarded as a first step in the study of dissipative 
processes in the universe utilizing the Full (non-truncated) theory. 

For a FRW (Friedmann-Robertson- Walker) cosmology the metric is given by 

dr 2 



ds 2 = -dt 2 + R(ty 



1 — kr' 



+ r 2 (d6 2 + sin 2 6d(f) 2 ) 



k = 0,±1, 



and the Einstein field equations and the energy conservation equation are given by 



H = -H 2 --(3 1 -2)p 
o 



-n, 



p = -3H( 1P + U), 
H 2 



1 k 



R 2 ' 



(1) 
(2) 
(3) 



where H = R/R is the Hubble expansion rate (we restrict ourselves to the expanding case 
only, i.e., H > 0), p is the energy density, and the local equilibrium pressure is assumed to 
obey 

P=( 7 -l)p, 1< 7 <2, 

with 7 constant. 

The bulk viscous pressure n obeys the evolution equation 



n = -3C# -rll- -rn 



3H 



t C T 
r~(~T 



(4) 



where ( > is the bulk viscosity coefficient, < r (= (j3 Q in |]]) is a relaxation coefficient 
for transient bulk viscous effects, and T > is the temperature. Equation (^) with e = I 
arises as the simplest way (linear in II) to satisfy the H-theorem (i.e., for entropy production 
to be non-negative ||). The Truncated theory effectively arises by setting e = 0, i.e., it 
corresponds to the case where the term in square brackets in equation (f|) is negligible in 



comparison with the other terms (see [T(J for the appropriate conditions). 

The Israel-Stewart theory is derived under the assumption that the thermodynamical 
state of the fluid is close to equilibrium, which means that the non-equilibrium bulk viscous 
pressure should be small when compared to the local equilibrium pressure, viz., 

|n|<p = ( 7 -i)p. (5) 

If this condition is violated, then one is effectively assuming that the linear theory holds 
also in the non-linear regime far from equilibrium. Such an assumption is unavoidable for 
viscous inflationary cosmology |8| . For a fluid description of the matter, equation (^) ought 
to be satisfied. However, note that non-linear viscous effects may arise in a phenomenological 



description of particle creation in the early universe ITT . 



II. THE DYNAMICAL SYSTEM 



Equations of state for ( and r and a temperature law for T are needed in order for the 



above system of equations to be closed. Belinskii et al. |T_2fl take ( and r to be proportional 
to powers of p, and this assumption is extended to T in 0. We shall follow [[jj and adopt 
'dimensionless' equations of state. That is, defining the dimensionless density parameter 

we shall assume that (/H and tH are proportional to powers of x, namely, 

jj = 3( x m , and T ~ W = hx n ) (7) 
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where m and n are constants which are assumed to be non-negative and ( Q and b are 
positive parameters. Clearly the equations of state employed will determine the qualitative 
properties of the models |]H,|5H3|,|12|| . Equations of state (0), which ensure that the asymptotic 



limit points represent self-similar models |T3|], are phenomological in nature and are no less 



appropriate than the equations of state used by Belinskii et al. []12||. We note that the 



equations of state chosen in [12 1 and those above coincide in the important case m = 1/2 = n 
(q = l/2 in!). 

From now on we shall take 

n — 0, a = b( Q . 

(Note that a, b are precisely the parameters used in [p]].) When n — 0, it follows that the 
relaxation rate is determined by the expansion rate: 

t- 1 = bH. (8) 

As argued in ||, for viscous expansion to be non-thermalising, we should have r _1 < H, for 
otherwise the basic interaction rate for viscous effects could be sufficiently rapid to restore 
equilibrium as the fluid expands. Therefore we impose the constraint 

b < 1 

on the relaxation parameter. 

Defining the dimensionless viscous pressure y and the new time variable t by 

n , dt TT 

y = — , and j t =H, (9) 

and using equation (P, equations (0) and (|J) become: 

x' =(x-l)[(3 1 -2)x + y], (10) 
y' = y[2 - b + y + (3 7 - 2)x\ - 9ax m - ~y*, (11) 

where 



and ' denotes a derivative with respect to t. Note that the linear condition @ becomes 

\y\ < 3(7 - l)x. 

Equations ( |I0D and ( [TT]) constitute a plane autonomous system of ODEs for x and y. In 
the Truncated theory e = 0, whence the final term in equation (O) is absent and there is 
no need to specify an equation for T. Hereafter we shall set e = 1, and adopt the following 
temperature power-law f|-f§: 

T = T // = T 3 r x r H 2r , with r = (13) 
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where the form of the exponent r follows from the integrability condition of the Gibbs 



equation when P = (7 — l)p [ IT0|, 15 1 - When the local equilibrium state of the expanding 
viscous fluid is thermalized radiation, then r = 1/4, in line with the standard Stefan- 
Boltzmann relation. Consequently 

H' x' 
* = 3 - 2(1 + r)— - (m + r)— , 
H x 

= Co + ciy + c 2 x + c 3 -, 

X 

where 

c = 5 + 2r + (37 - 2)(m + r), 
ci = 1 — m, 
c 2 = (37 - 2)ci, 
C3 = m + r. 



A. Flat universe 

All of the FRW models are governed by equations ( [XQ|) and (|11]) together with equation 
(0). We note from (|T0D that 2 = 1 is an invariant set, where from (||) we see that this set 
represents the flat FRW models. Let us study this physically important zero-curvature case 
first. When x — 1, the thermodynamic laws are simplified to (§), (fl3|) and, by (0), to 

5 



C « H. (14) 

Thus the bulk viscosity coefficient, like the relaxation rate, is also determined by the expan- 
sion rate. Furthermore, 

* = (co + ca) + (ci + c 3 )y, (15) 



whence equation becomes 



y' = -^—^y 2 -by-9a. (16) 

That is, the equations governing the evolution of the flat FRW viscous fluid models reduce 
to a single autonomous ODE in y. Since < r < 1/2, ( PS|) is a Riccati equation with 
constant coefficients and its solutions can be found in implicit form. 
Defining the positive parameter 

B x = b 2 + 18a(l -r), 

it follows that there are two equilibrium points, one positive, one negative, 

y ± = 1 (b±JF 1 ), (17) 

where one is a sink and the other is a source (with respect to the invariant set x = 1, not 
the full set of all FRW models). The points (fL7f) correspond to the special solutions found 



in 



16] and re-discovered in pj 
Therefore, the behaviour of the flat models using the Full (non-truncated) theory is qual- 
itatively the same as the behaviour in the Truncated theory |l[]. Of course, this qualitative 
similarity only holds for the restrictive thermodynamic laws (B), (|T3|) and fll4|). 



B. Curved universes 



Let us now return to the general curvature case i ^ 1 [see equations ( |IDD and (|TT|)]. 
Equation (|TTD can be written as 
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2 + f) + *(*-2)(f-l) + ^-l) + ^ 



-1 



9ax m . (18) 



There are two equilibrium points lying in the invariant set a; = 1, namely (1, y^) where 
y ± is given by equation fllTf ). Previously (in the case of the flat models) we considered the 
stability of the equilibrium points only with respect to the invariant set x = 1; let us now 
discuss the stability of these equilibrium points with respect to the curved FRW models. 
The equilibrium point (l,y + ) is a source with the invariant set x — 1 as one of its primary 
eigendirections. If y~ + 37 > 2 then the equilibrium point (1, y~) is a saddle with x = 1 as 
the stable manifold. If y~ + 37 < 2 then the equilibrium point (1, y~~) is a sink and hence it 
represents a future asymptotic attractor. From equation (|IDD, the equilibrium points (x,y) 
not lying in the invariant set x = 1 satisfy y = —(37 — 2)x, and hence from flTBD we obtain 



9ax m -^(37-2)(26 + 2r + l)x = 0. (19) 

For m > 0, there exists a singular point at the origin (0, 0). [Note, however, that the system 
of ODEs as given by equations (|10D and (|1^) is not defined at x = except when c 3 = 
and therefore the point (0, 0) may not be a well defined equilibrium point of the system.] 
Changing to polar coordinates, it can be shown that this singular point is saddle-like in 
nature (hyperbolic sectors) if m < 1. If m > 1, then the point (0,0) has parabolic and 
hyperbolic sectors. 

If m ^ 1, then there is a second equilibrium point at 

l/(m-l) \ 



If 



-(37-2)5 . (20) 



fi 2 = (37 - 2)(26 + 2r + 1) - 18a > 0, 

then x > 1, and when m < 1, this point is a saddle. If i?2 < 0, then 5 < 1, and when m > 1, 
this equilibrium point is again a saddle. There is a variety of other possible behaviours. 



III. DISCUSSION 
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A. Exact Solutions and Asymptotic Behaviours 



The qualitative behaviour of the flat FRW models has been determined completely. The 
unphysical flat models evolve from the equilibrium point y = y + at i = — oo, where y = y + 
corresponds to the solution (after recoordination) 

R(t) = R t^y ++ ^\ H(t) = ^—t-\ 

y+ + 37 

towards either points at infinity or to the point y = y~ (at i = —00), which, if y~ 7^ —37, 
has solution 

R(t) = Ro(t - t ) 2/( ^ +37) , H(t) = —^-(t - U)- 1 

y- + 37 



[Note that if y + 37 > 0, then the solution (|22|) can be recoordinatized such that t = 0. 



These models and the equilibrium point y = y + are unphysical since they have positive bulk 
viscous pressure. Those models which evolve towards y = y~ have negative bulk viscous 
pressure after a certain time, and may be considered as physical models after this time. The 
models which are physical for all times (i.e. which have II < for all times) are (^) and 
those which evolve from infinity at t = —00 towards y — y~ at i — 00." [Note that, by @, 
(pip and (|22D, t = — 00 corresponds to t = 0, while i = 00 corresponds to t = 00.] 



If y = —37 then the solution has the form 

R(t) = R e Hot , H{t) = H 

P (t) = ml II(t) = y -H, 2 . (23) 

The exponential inflationary solution (^) clearly violates the condition (|5|) (cf. 0). The 
solution ( |22"D violates ([5]) if y~ + ?> r y < 3, when the expansion is driven by a large and effective 
nonlinear bulk viscous pressure. The expansion is from a big bang, and is more rapid than 
in the corresponding equilibrium solution (y~ = 0). Indeed, if 
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y~ + 3 7 < 2, (24) 

then the solution represents a power-law inflationary solution. Knowing that condition (|2"1D 
is also the requirement that the equilibrium point (1, y~) be stable, we can conclude that the 
power-law inflationary solution (^) is the future asymptotic attractor for all bulk-viscous 
inflationary FRW models. 

We emphasise that bulk viscous inflationary solutions, such as fl2"2|) , violate the condition 
(H), so that their existence is dependent on assuming that the theory holds in the nonlinear 
regime. Furthermore, these inflationary solutions are limited by the simple equation of state 
P = (7 — so that, in particular, they cannot account for the processes necessary to 
provide an exit from inflation. The solutions are at most valid during inflation, and more 
realistic models would be needed to incorporate exit and re-heating. 



B. Entropy Production 

On physical grounds, one expects that y < 0, since the evolution of specific entropy is 
given by §,§ 

s = — (25) 

where n is the number density. We note that solution (|22"D always satisfies y < 0. From 
equation fl25|), it follows that the growth of entropy in a comoving volume between times 
to < t < t\ is given by 

- Efo) = ~ ^dt, (26) 

where k is the Boltzmann constant. The amount of entropy generated can be calculated for 
each of the solutions (^l|) , (|22|) , and (^) . Analyzing the physically more relevant case (^2|) . 
we find that (re- instating constants previously set to unity), 

YU \ YU \ - ^ 3(1 ~ r)c2 ( R 3 \ ( 2 V 1 /, -%-/[7(2/"+37)] + -2»-/[7(«-+37)] 



(27) 



where c is the speed of light and G is the gravitational constant. By equation (|22"D , we must 
have y~ + 37 > for an expanding solution. This ensures that ( p7|) gives X(ti) > S(t ). 

The amount of entropy generated in the de Sitter model (solution (|23|)) is the same as 
that calculated by Maartens jSJ, namely 

- E(t„) = ( #1 #0^ (e 3Hotl - e 3Hoto ) . (28) 



8ttGA; V T / 

It is shown that bulk viscous inflation can generate significant amounts of entropy without 
re-heating. For the de Sitter model, using typical parameters of inflation, and assuming 
that almost all of the entropy is produced by inflation, one finds the following value for the 
amount of entropy produced during exponential inflation || 



S^2.1xl0 87 , (29) 



which is in agreement with the expected value. The power-law inflationary solution ( |22| ) 
[with y~ + 37 < 2, i.e., satisfying equation (p4"D1 has less efficient entropy production, but 
nonetheless can also produce significant amounts of entropy. 

In the above, we have only considered entropy production in the models corresponding 
to the equilibrium points of the dynamical system. By considering a simple example, we 
can investigate the entropy production in the more general flat FRW models. We choose 
parameter values r = 1/4 (necessarily 7 = 4/3), a = 1/27, and 6=1/4. In this case the 
differential equation (|16D reduces to 



, 3 

V = a 



(30) 



which has a solution of the form (neglecting the constants of integration) 

J I - tanh(|f) \y - || < 1 
V = I I - coth(ft) |3/ — || > 1 

The Hubble parameter is given by 



(31) 



e"~ r cosh 4/3 (|t) |y-||<l 

H=\ * t - • ^4/3,3* L. 1, . 1 ( 32 ) 



c -' smh 4/3 (|t) \y- || > 1 
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and R = Roe*. The change in entropy in a comoving volume produced between times 
to < t < i\ is then given by 



Eft) - E(t ) 



±2(e- fl/4 - e~ h/A ) + (e"* 1 - e"'" ) + (e 



fi/2 _ to/2 



e t0/2 ) . (33) 



v 167tGA;To 

It can be concluded in this simple model for ii/io > 1 that as i\ increases the entropy in a 
comoving volume grows exponentially with respect to t\. 



IV. CONCLUSIONS 

The behaviour of the viscous fluid FRW models where the bulk viscous pressure satisfies 
the Full Israel-Stewart theory of irreversible thermodynamics has been analyzed. The sta- 
bility of the equilibrium point (0, 0) representing the Milne model depends upon the value 
of m which appears in the equation of state for the bulk viscosity. The equilibrium point 
(x, y) can represent either an open, flat or closed FRW model depending upon the value of 
the parameter B2. Exact determination of the nature of this particular equilibrium point 
is extremely difficult. However, a partial result is possible: if -82(1 — m) > 0, then the 
equilibrium point is a saddle. There exist two equilibrium points with qualitative behaviour 
similar to that found using the Truncated Israel- Stewart theory. 

It can be concluded that the behaviour of the FRW models in which the bulk viscous 
pressure satisfies the Full Israel- Stewart theory can in principle be qualitatively similar to 
the behaviour of the FRW models in the Truncated theory. One cannot say, however, that 
the Full theory has the same behaviour as the Truncated theory in all cases because it is not 
at all clear what effects the presence of anisotropies or different equations of state will have. 
For example, in the models studied here, it was the equations of state for the temperature 
and for the bulk viscosity coefficient that played major roles in determining the dynamics 
of the models. In the case of a relativistic Maxwell-Boltzmann gas, which has very different 
equations of state, the Truncated and Full theories can lead to very different behavior, with 
the Truncated theory leading to pathological behavior of the temperature in many cases . 
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As stated above, the consistency condition that viscous expansion should be non- 
thermalising requires b < 1. Further constraints may arise from entropy arguments. The 
evolution equation (f|) already guarantees that entropy production is non-negative. But 
one may place constraints on the rate and amount of entropy production. If we impose the 
requirement that the specific entropy, s, should increase with expansion, but at a decreasing 
rate, then we have y < and possibly further constraints on r (equivalently 7) and m. 

A complete analysis of the asymptotic behaviours of these viscous fluid models depending 
on the (many) free parameters in the model (a, b, 7, m) and utilizing the energy conditions 
can be made. However, the next step in this research programme is to attempt to use results 
from kinetic theory in order to motivate physically plausible equations of state, or, at the 
very least, to limit the form of the phenomological equations of state used. 
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